A Magnetic Resonance Realization of Decoherence-Free Quantum Computation 
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We report the realization, using nuclear magnetic resonance techniques, of the first quantum 
computer that reliably executes an algorithm in the presence of strong decoherence. The computer 
is based on a quantum error avoidance code that protects against a class of multiple-qubit errors. 
The code stores two decoherence- free logical qubits in four noisy physical qubits. The computer 
successfully executes Grover's search algorithm in the presence of arbitrarily strong engineered 
decoherence. A control computer with no decoherence protection consistently fails under the same 
conditions. 
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A computer that uses the laws of quantum mechan- 
ics to store and manipulate information could in theory 
perform certain tasks such as searching 1] and factor- 
ing 2] with incredible efficiency. The most critical prob- 
lem that must be solved to make quantum computing 
possible on a useful scale is decoherence, the inevitable 
process of entanglement between a quantum computer 
and its environment. Decoherence causes the superposi- 
tion states that carry information within the computer 
to decay rapidly. Several solutions to the decoherence 
problem have been proposed (for a review, see Q). One 
technique, quantum error avoidance, calls for information 
within the computer to be carried exclusively by quan- 
tum states that are not adversely affected by decoher- 
ence 0,03, 00 (f° r a review, see @). Here we present the 
first experimental proof that a nontrivial quantum com- 
putation can be protected against decoherence Us- 
ing quantum error avoidance, we have constructed a nu- 
clear magnetic resonance quantum computer 0, 0, 
which is unaffected by certain types of decoherence. Our 
computer successfully executes Grover's quantum search 
algorithm [jj in the presence of arbitrarily strong engi- 
neered decoherence. A control computer with no deco- 
herence protection consistently fails under the same con- 
ditions. 

Decoherence is typically characterized by the decay 
of off-diagonal elements in a system's density matrix p. 
Formally, decoherence takes the system from state pi 
to a state pf = ^2dEdPiE\ where the Kraus operators 
Ed describe transformations that may result from the 
system-environment coupling (they satisfy ^ d E^Ed = 
7 = identity) HE1- 

When the coupling between a quantum system and 
its environment possesses an element of symmetry, some 
of the system's states will be immune to decoher- 
ence 0, IE IS 13 • These states span a decoherence- free 
subspace (DFS). The quantum computer we have con- 
structed comprises two decoherence-free logical quantum 
bits (qubits encoded in the DFSs of four noisy phys- 



ical qubits. The code protects against multiple-qubit 
errors To satisfy the symmetry condition for the 
existence of DFSs, we assume the system-environment 
coupling affects certain pairs of qubits rather than af- 
fecting each qubit independently. Note that this error 
model is different from the popular "collective decoher- 
ence" model HUB]. 

The multiple qubit errors model is relevant to a number 
of physical systems recently used as quantum comput- 
ers [3. For example, the primary source of decoherence 
in liquid state NMR is the random modulation of internu- 
clear dipolar interactions by molecular tumbling. Under 
certain conditions, for example very slow tumbling, the 
interactions reduce to a symmetrical multiple qubit er- 
ror process. This makes some of a system's coherences 
resistant or immune to decoherence (an effect recently ex- 
ploited in NMR studies of large proteins |lj| ) and gives 
rise to DFSs similar to those used in this work. However, 
the object of this work is to demonstrate DFS protection 
against multiple qubit errors, not to demonstrate specific 
resistance to the natural decoherence processes of liquid 
state NMR. We chose an error model that supports a rel- 
atively simple DFS and that affords us complete control 
of the decoherence strength, and as such it is not related 
to our system's natural decoherence. 

Specifically, the code our computer uses resists errors 
of the form [15|. 

Ed — ddfilihhh + a<j, 1-^1X2/374 

+ ad^IihX^Xi + ad,3X 1 X 2 X 3 X4, (1) 

where X n indicates that physical qubit n € {1, 2, 3, 4} is 
flipped and I n indicates that it is unaffected. There are 
four DFSs for this set of errors |j| • Each DFS is a simulta- 
neous eigenspace of the operators {I1I2I3I4, X1X2I3I4, 
/1/2A3A4, X 1 X 2 X 3 X i } with eigenvalues ±1. The fol- 
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lowing states are an orthonormal basis for one DFS: 

|00)i = (|0000) + |1100) + |0011) + |llll))/2 

|01)i = (1 1000) + |0100) + |1011) + |0111)) /2 

|10)i = (|0001) + 1 1101) + |0010) + |H10))/2 

= (|1001) + |0101) + |1010) + |0110))/2. 

These can be used as basis states for two decoherence- 
free logical qubits and are labeled as such. The other 
three DFSs are related to this DFS by sign changes. For 
example: 

|00)| = (|0000) + |1100) - |0011) - |llll))/2 
|00)| = (|0000) — |1100) + |0011) — jllll)) /2 
|00)| = (|0000) - |1100) - |0011) + |llll))/2, 

where superscript numerals indicate to which DFS a state 
belongs. The remaining basis states of DFSs 2-4 are 
obtained by applying similar sign changes to the other 
states of DFS 1. 

The code uses all four of these DFSs in classical paral- 
lel. An arbitrary state of the two logical qubits, — 
a|00) L + 6|01) L + c|10) L + d|ll) L where |a| 2 + \b\ 2 + |c| 2 + 
\d\ 2 = 1, is encoded in the density matrix p\^) L describing 
the four physical qubits according to 
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p mL Htf>iMz=I>MiMi, (2) 

i=l 

where \ip)i = a|00) J L + 6|01)^ + c|10) l L + d|ll) l L and 
ci = 1/4. Note that only the component of the den- 
sity matrix that deviates from identity is described above 
(the identity portion of p is immutable and unobservable 
during any NMR experiment 0). 

It should be noted that in general a quantum superpo- 
sition of states from different DFSs is not decoherence- 
free. This is because the DFSs are eigenspaces of the 
Ed operators, and for a given E d each DFS may have a 
different eigenvalue. However, the encoding of Eq. 
uses the DFSs in classical superposition only and P\$) L is 
therefore unaffected by the errors described by Eq. QJ. 
It can be shown that E^iPYl = ^£mI'</')l where Ad.i is a 
constant @- It follows that, when the errors of Eq. (JJJ 
occur: 
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Pf = Y, E *p E l=Y,Y, c * E ^)^iE d 

d d i— 1 

i=l V d ) i=l 

where c! i = c, J^d Ad,iA d i . Thus the net effect of the er- 
rors is a change in the relative weights c[ of the different 
DFSs: The logical qubit information in each DFS is in- 
tact. For the errors we implement experimentally, it is 
always the case that ^2 d A d ^A d i = 1 so that the errors 



have no effect whatsoever {d i — Cj so that pf = p). This 
simplifies interpretation of the experimental results but 
is not essential to the code's performance. 

Pulse sequences that perform logic gates on the two 
encoded qubits were developed using methods derived 
in 0| and will be described in detail in a subsequent 
publication. Unfortunately the computer leaves the DFS 
code during gate sequences. It has been shown that this 
class of DFS codes can function as active quantum er- 
ror correction codes against errors occuring during qubit 
manipulation [l^ | , a property that can be used in future 
implementations to detect and correct such errors. 

Before any computation, the computer's two logical 
qubits are initialized by temporal averaging [l7| to the 
state |00)j> The code maps this logical qubit state to 
the following state of the four physical qubits: 

P\oo) L = (|0000) (0000|+ |1100)(1100| 

+ |0011)(0011| + |llll)(llll|)/4 (3) 

This density matrix can be decomposed into a sum of 
tensor products of the 2x2 identity matrix / and the 
Pauli matrix Z . 

P\00) L — (Ilhhl4 + Z1Z2I3I4 

+ hI 2 Z 3 Z 4 + Z 1 Z 2 Z 3 Z 4 )/16 

The I1I2I3I4 term can be neglected and each of the 
remaining three can easily be prepared from the sys- 
tem's e quil ibrium state using standard pulsed NMR tech- 
niques 18] . During a quantum computing experiment, 
the computation is repeated three times, each time pref- 
aced with a pulse sequence that prepares a different one 
oiZxZ^I-j,!^ I1I2Z3Z4, and ZiZ 2 Z 3 Z 4 . The three results 
are added together and because computation is a linear 
quantum operation, the summed result corresponds to a 
computation starting from the state |00)l. 

Because the code uses states from four DFSs in clas- 
sical parallel, our ensemble computer does not use true 
pseudo-pure states. We chose this approach to reduce the 
number of temporal averaging steps, facilitating a thor- 
ough test of the computer's resistance to decoherence. In 
fact, the three temporal averaging steps we perform are 
a subset of the fifteen required for a pseudo-pure state 
implementation. 

We use our computer to perform Grover's quantum 
search algorithm [l| (Fig. QJ. The algorithm's purpose 
is to retrieve, from an unsorted list, the single item 
that satisfies a given criterion. Grover's algorithm is 
highly efficient, requiring only 0(y/~N) steps to search 
a list with N items; a classical algorithm requires O (N) 
steps. In our implementation of Grover's algorithm, the 
logical qubit basis states |00)l, |01)l, |10)l and |11)l 
correspond to the four items in our list. We choose 
state 1 11) z, to correspond to the item we wish to re- 
trieve. The algorithm's first step is to prepare the regis- 
ter in an equal superposition of its basis states, \iP)l = 
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FIG. 1: Running Grover's search algorithm on an error-avoiding quantum computer. Time runs from left to right. Decoherence- 
free logical qubits are represented by horizontal lines. Engineered decoherence is applied at nine points during the experiment 
to test the computer's resistance. 




FIG. 2: Isotope-substituted glycine molecule. Spin- 1/2 nuclei 
used as physical qubits are numbered 1-4. The chemical shifts 
of 13 C spins 1 and 3 differ by 16.5 kHz on a 500 MHz ( X H 
resonance frequency) NMR spectrometer. 

(V 2 ) X^ 1= u X^ 2 =o \ x i x z)l- The rest of the algorithm is 
an iterative process that increases the amplitude of the 
sought-after state (|11)l) until it is the dominant part of 
the superposition. For our two-qubit register, only one 
iteration is required. Finally, the logical qubits are mea- 
sured; an outcome of indicates that the algorithm 
has been successful. (We have also used our computer 
to implement the improved Deutsch-Jozsa algorithm de- 
scribed by Collins et al. 0], with similar results.) 

Our computer is based on the isotope-substituted 
glycine molecule shown in Fig. [21 20]. Each of the 
molecule's four spin-1/2 nuclei is used as a physical 
qubit. The errors of Eq. are applied artificially ac- 
cording to the following protocol. At each of nine points 
in the experiment (indicated in Fig. Q, error operator 
X1X2I3I4 is applied with probability e, then error op- 
erator I\ I2X3X4 is applied with the same probability. 
To simulate the effect of a microscopic random process, 
the experiment is performed 2048 times, each time with 
different randomly selected errors, and the results are 
averaged, giving the overall decoherence process a non- 
unitary, deterministic character. The resulting decoher- 
ence increases with the error probability e, becoming 
strongest at e = 0.5. Formally, the operators describing 
our engineered decoherence are Eo — (1 — e) ii I2I3I4, 
E x = ^e(l-e)X 1 X 2 I 3 h, E 2 = (1 - e)lil 2 *3*4, 
and E3 — 6X1X2X3X4. (Note that for these error op- 
erators, it is clear that ^2 d EdP\m) L E\ — p\m) L , where 



p\oo) L is defined by Eq. J3J).) 

We have repeated the Grover algorithm experiment in 
the presence of nine different levels of engineered deco- 
herence ranging from e = to e = 0.5. At all values 
of e, the resulting NMR spectra contain little distortion 
and clearly describe the final state of the qubit register 
as indicating the computer has successfully exe- 

cuted the algorithm. To quantify the resistance of each 
computation to the applied decoherence, we measured 
the integrated absolute intensity of the final signal rela- 
tive to the e = experiment. The dependence of signal 
intensity on e is different for each of the experiment's 
three temporal averaging steps and we have chosen to 
analyze the results of the steps separately (Fig. EJ). Wc 
observe only small losses of signal, and these cannot be 
attributed to any fault in the DFS encoding. The losses 
are predominantly due to imperfections in pulses used to 
implement the engineered decoherence, evidenced by a 
linear decrease in signal with increasing e. 

As a control, we have repeated these experiments on 
a quantum computer that docs not use error avoid- 
ance. The unprotected computer is similar to the error- 
avoiding computer, with only the following changes: 
Spins 1 and 4 (Fig. [2Jl serve directly as qubits (in place 
of the two logical qubits), the temporal averaging scheme 
prepares a different (pseudo-pure) initial state, and dif- 
ferent pulse sequences are used to implement quantum 
logic gates. The same engineered decoherence is applied 
and the computer executes the same algorithm as in the 
error-avoiding computing experiments. As expected, the 
unprotected computer's signal intensity decreases rapidly 
with the error strength e. We observe effectively no signal 
when e > 0.3 and the result of the algorithm is incorrect 
or unreadable for e > 0.2. 

The DFS encoding our error-avoiding computer uses 
has an overhead cost, but the results show it is small 
compared to the protection it affords. The pulse se- 
quences that perform logical operations on the error- 
avoiding computer are more complex than for the unpro- 
tected control computer, so the error-avoiding computer 
is more vulnerable to signal loss due to pulse imperfec- 
tions and natural spin relaxation. This is why, when sig- 
nal intensities from the two computers are compared on 
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FIG. 3: Experimental results in the presence of decoherence, for both the error- avoiding and unprotected computers running 
Grover's algorithm. Results from the three temporal averaging steps appear separately in (a, Z1Z2I3I4,), (b, Z1Z2Z3Z4), and 
(c, I1I2Z3Z4). The integrated absolute intensity of the output signal measured at the end of Grover's algorithm is charted as a 
function of decoherence strength e. For some values of e, the signal from the unprotected computer is phase inverted with respect 
to the correct output signal: In such cases we report a negative signal intensity. Square and round points are experimental 
data from the error-avoiding and unprotected computers, respectively, with each dataset normalized to the intensity of its first 
point. Solid lines are linear regressions. Dashed curves are a theoretically predicted intensity function, S = (1 — 2e) n with 
n = 6 in (a, c) and n = 12 in (b), obtained by predicting the unprotected computer's state at each of the experiment's nine 
decoherence points and counting the number n of error operators which could change the system's state. Dotted lines are signal 
intensities from the error-avoiding computer normalized to the intensity of the unprotected computer's e = signal; the signals 
in (a) were recorded on different nuclei and cannot be compared in this way. 



an absolute scale (dotted and dashed lines in Fig. |3J) , the 
unprotected computer gives the stronger signal for low 
values of e. However, the unprotected computer consis- 
tently fails when e > 0.2, while the error-avoiding com- 
puter gives the correct result for all e. The overall fidelity 
of the Grover algorithm is governed by the temporal aver- 
aging step that is least tolerant to decoherence (Fig. EJ3) , 
and for this step the error-avoiding computer's signal is 
the more intense for e > ^0.15. Even at this low level 
of decoherence, the protection afforded by DFS encoding 
outweighs the overhead involved. 

In summary, we have provided the first experimental 
demonstration of quantum computation in the presence 
of strong decoherence , thus proving that quantum er- 
ror avoidance based on DFSs can very effectively protect 
qubits from decoherence during the execution of a quan- 
tum algorithm. We have implemented Grover's search 
algorithm on two two-qubit quantum computers, one 
error-avoiding and one unprotected. While the unpro- 
tected computer fails when exposed to even a moderate 
amount of decoherence, the error-avoiding computer is 
successful in the presence of the strongest possible deco- 
herence. This demonstration is a proof of the concept 
of quantum error avoidance and suggests that DFS en- 
coding will play an important role in future experimental 
implementations of quantum algorithms in the presence 
of decoherence. 
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